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Abstract 

The thermodynamics of interacting vortices constricted to move in one-dimensional 
tracks carved out on superconducting films in the extreme type-II limit is mapped 
into the Coulomb gas model of random matrices. Application of the Selberg Inte- 
gral supply exact expressions for the constitutive relation for two configurations of 
tracks. 
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In the present paper we derive an exact expression for the constitutive relation be- 
tween the external magnetic field, H, and the magnetic induction, B, inside a type II 
superconductor, i.e., we determine the relationship B = B(H). The induction inside the 
superconductor is determined by the density, n, of magnetic vortices. Each vortex carries 
one flux quantum <& , hence B = $o n -H m order to be able to find an exact expression 
for B(H) we restrict ourselves to thin superconducting films. The interaction energy of 
two vortices depends in this case logarithmically on the distance of separation of the two 
vortices. Furthermore, we assume that the vortices are confined to move on tracks en- 
graved into the superconducting film, see Fig. 1. This experimental situation has in fact 
been studied by Pruymboom et al.pf The thickness of the film is reduced in the track. 
Since there is a loss of condensation energy in the core of a vortex, it will be energetically 
favourable to position a vortex with the core on the track rather than outside the track, 
where the core would have to drill its way through a thicker layer of superfluid. 

The combined simplification of logarithmic interaction and restriction on the config- 
urational degrees of freedom allow us to map the problem onto the Coulomb gas model 
of random matrices. The partition function is known exactly for this model and the 
constitutive relation readily derived. 

Let us now describe the derivation of the relation between H and B. The partition 
function for a fixed number of N two dimensional vortices on a track is given by: 

t N r L dx \ 

Z N = [R I exp (-0E[*i, • • , (1) 

where 

E[x 1 ,...,x N ]=NE e -E 1 £ \n( lXa ~ Xbl ), (2) 

where E c (> 0) is the energy of the vortex core and E\(> 0) sets the scale of the vortex 
interaction. Both E c and E% are temperature dependent. 0] The length of a the tracks 
are denoted by L, see Fig. 1. The scale of the logarithmic interaction is set by the 
effective magnetic penetration depth, A e ff, of the film. Expressed in terms of the bulk 
penetration depth A and the film thickness d one has A e ff = X 2 /2d. p] (3 = 1/T is the 
inverse temperature. N is related to the magnetic induction B through BLR = N& , 
where R is the distance between the individual tracks, see Fig. 1. / the length scale of the 
configuration space and may be taken to be equal to £(0), the coherence length at zero 
temperature. We assume that L < X c s and R 3> A e g. 
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The condition A c g ^> L ensures that the logarithmic interaction is valid for all pos- 
sible separations of a pair of vortices on the same track. Beyond the distance A e g the 
vortex interaction decays exponentially due to the diamagnetic screening currents. [|], [5] 
Consequently, we can neglect the interaction between vortices on different tracks if the 
condition R ^> A e g- is fulfilled. 

It should be noticed that the vortices will stay in the tracks for magnetic fields less 
that a certain filed strength H*. At H* the repulsion between the flux lines becomes so 
strong that it will be able to over come the pinning force exerted by the carved track. The 
actual value of H* will depend on the depth of the groove and in a first approximation 
will be equal to the H Cl of the bulk material. The partition function in Eq. 1 applies for 
fields H < H*. 

The partition function can be re-written as: 

- /L\ N f L \ N ( N -i)Ei ( n fX \ 
^v = e"Wy) (f-) [U dt a ) n \ta-h\ 2E \ (3) 

where we have introduced the following dimensionless energies: 

E = — E = — (A) 

and the dimensionless integration variables t a . The multiple integral in Eq. [3] is a special 
case of the Selberg integral, with its aid we obtain the partition function and the free 
energy of this system: 

F N 



.-^=\nZ N = N 



\n(L/l) - lnr(l + E x ) - E c \ + iV(iV - l)£iln(L/A eff ) 

x rr(i + 0' + i)A)r 2 (i+^0" 



j=0 



(5) 



T(2 + (N-l+j)E 1 ) 

An expression for the lower critical field H Cl is obtained from the condition that the 
Gibbs free energy of the state with a single vortex F± —H C1 BV/An (V = LRd is the volume 
per track) is equal to the Gibbs free energy of the state without any vortices F = 0.0 
The result is 

H cl = ^(E c -Tln(L/l)) (6) 

The constitutive relation is obtained in the following way which makes use of the fact 
that H acts as the chemical potential for the vortices and ^ = 0, where G is the Gibbs 



3 



free energy for the system: |7]] 

dF N /T _ 1 $o H 
dN ~ A^TrY 



V (7) 



We simplify the above expression in the large N limit, by replacing the sum of the log- 
arithms by an integral while discarding the differences between N and N ± 1. With the 
aid of T(ax + b) « (27r) 1 / 2 e - as (aa;) a:r+f, " 1/2 , a > 0, 6 > 0, ac > 1 we find 

^U = 2X ln(4A eff /L) + lnX, (8) 

47TJ 

where X := NE\ = (BLR/<& )Ei, and we have discarded the N independent and 0(A^ _1 ) 
terms. The magnetic permeability fi in the high field limit, can be found by further 
discarding the lnX term 

B = HoH, (9) 

with 

47rE 1 L J Rln(4A eff /L)' (10) 
The constitutive relation has a weak logarithmic correction for not too large H. This is 
found by the substitution X = X Q + SX, where 

/ dH$ \ 

6X i / 8vrTln(4A cft /L) , v 

^ I l + f^ J' 1 ' 

which in turns supplies a non-linear permeability n(H) = fi (l+SX/X ). The permeability 
is slightly reduced from the its high field values which is to be expected. 

So far we have considered the ensemble where the position of the vortices are restricted 
to a straight line of length L. Let us now consider a system in which the vortices are 
confined to a circular track of radius R Q . In this case we avoid edge effects. The ther- 
modynamics of the the logarithmically interacting vortices on a ring are obtained from 
Dyson's circular ensemble. |]] The interaction energy of the vortices is given by 

E C [6 1 ...9 N ] = NE C -E 1 ln(^|e^-e^|V (12) 

A - 



l<a<b<N 

and the partition function 

N 



cS 



( r 2n R \ 

Z n={J[J o ^r l d6 a jexp(-f3E c [6 l ...e N }), (13) 
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where R (<^ A e g) is the circumference of the circle. The constitutive relation is in this 
case easily obtained directly from the difference F N+1 — F N in free energy. We have 

( F N+1 - F n\ 



T((N + 1)E 1 



= 2NE 1 \n(R /\ eS )+\n > . (14) 

For large N and with the aid of large x behaviour of T(ax + 6), we find, 

®" H V = 2\n( y \ cS /2nR )Y + E^nY, (15) 



87r 2 R dT 
where 

Y:=NE U (16) 

The result for the permeability in the high field limit reads 

1 $ 2 

n c = Z ^ (1*7) 

^° 327r 2 E 1 R ]n(\ e s/2irR ) 1 ; 

a result similar to that of the straight track. 

We have derive an exact expression for the constitutive relation for an artificially pat- 
terned superconducting thin film. The result should be directly experimentally accessible. 
Pinning along the the engraved track may in principle influence the motion of the vortices. 
Although, it appears that for deep grooves pinning should effect the motion parallel to 
the track rather insignificantly. 

Our future extension of the present calculation will be to apply random matrix theory 
to the dynamics of the vortices in the constricted geometry. 
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Captions 



Figure 1. 

Sketch of a superconducting thin film with tracks engraved into the surface. 



Figure 2. 

The induction B as function of the external magnetic field H obtained from Eq. 8 for 
the following set of parameters: L = 1/xm = R/10, d — 10 A, and A = 1500 A. The 
temperature is supposed to be 10K. 
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